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^j i Abstract 

Muon spin rotation/relaxation/resonance (MuSR) technique for studying matter structures is consid- 
ered by means of a recently introduced probability representation of quantum spin states. A relation 
,-G between experimental MuSR histograms and muon spin tomograms is established. Time evolution of 

muonium, anomalous muonium, and a muonium-like system is studied in the tomographic representa- 
tion. Entanglement phenomenon of a bipartite muon-electron system is investigated via tomographic 
analogues of Bell number and positive partial transpose (PPT) criterion. Reconstruction of the muon- 
electron spin state as well as the total spin tomography of composed system is discussed. 

cr: 
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Introduction 



The description of spin states is usually based on the concept of density matrix. The density matrix 
formalism was introduced originally by Landau pQ and von Neumann [2J for continuous variables, and 
then successfully extended to the case of discrete variables and spin states (see, e.g. [3j). In contrast to 
the state of classical particle, the density matrix cannot be measured immediately in the experiment, 
as the outcomes of quantum observables have a probabilistic nature. The attempt to bring classical 
and quantum descriptions closer together was made by Wigner [3] , who introduced the quasi-probability 
function W(q,p) in phase space such that marginal distributions J W(q,p)dp and J W(q,p)dq are true 
probability distribution of position and momentum, respectively. However, Wigner function can take 
negative values. Moreover, there is no direct way to measure VF-function. So far many quasi-probability 
functions were introduced, for instance, Hisimi Q-function [5] and Sudarshan-Glauber P-function [BJ 
[7J (the review of different quasi-probability distributions is presented in [8]). Experimental needs to 
reconstruct density matrix or Wigner function of photon states by means of experimentally measured 
probability distributions resulted in developing tomographic methods (both theoretical [9llll][l~8] and 
experimental [12013] ) based on Radon transform |14pi5j . Though tomograms were originally considered as 
auxiliary tools, recently it was proposed to treat tomograms as primary notions of quantum states |16pi7j 
because a tomogram is a fair probability distribution function, contains the same information about the 
system as a density matrix, and the crucial fact is that tomograms can be measured in laboratory. The 
tomographic description of light states was extended to include optical tomography 0CEO], symplectic 
tomography [IT], photon-number tomography [18-20], and their recent generalizations [21[l22j. The 
tomographic-probability representation was then extended to spin states [23,24] and developed in [25][26]. 
In the paper [27], the general form of unitary spin tomography was considered. The spin tomography 



with a finite number of rotations was discussed in [28J . We refer the reader to the review [29] for detailed 
history of the tomographic-probability representation. 

The aim of this paper is to apply spin tomograms to muon spin rotation/relaxation/resonance ex- 
periments usually abbreviated to MuSR or acronym /xSR. At present, MuSR method is primarily used 
to get insight about the microscopic behavior of materials and nanostructures (see the reviews |30H32j ) . 
although the origin of this method is due to experimental verification of the proposal that the weak 
interaction might violate parity symmetry, and measuring the anisotropy of muon decay predicted by 
V-A theory of weak interactions. Muons (//) are utilized in //SR technique as probes. Muon is a short- 
lived subatomic particle and has the following properties: charge ±1 elementary charge (fi + and fj,~ , 
respectively), mass m^ = 206.77 m e , spin s = 1/2, magnetic moment /i M = 3.18334 fi p , mean lifetime 
r = 2.197 • 10~ 6 s. Positive (negative) muon decays into positron (electron) and a neutrino-antineutrino 
pair as follows: 

/i + — >• e + + u e + Ufj,, \T — >e~+i) e + v tl . (1) 

Thanks to violation of parity symmetry, during the decay of /i + , the positron e + is preferentially emitted 
along the direction of muon spin. The decay positron can be emitted along an arbitrary direction but 
the probability of such an event depends on the direction n = (cos <p sin 8, sin if sin 9, cos 0) (anisotropy). 
Nevertheless, it is possible to measure the angular distribution T(n) of the decay positrons from a bunch 
of muons (ensemble) deposited at the same conditions. If detectors are capable to detect positrons 
of an arbitrary energy (from to oo) and muon spin is aligned along the z-axis, i.e. its vector state 
|j(w = 1/2, m^' = 1/2) is an eigenstate of angular momentum operators J z and J 2 , then the angular 
distribution T(n) reads (see, e.g. [30J) 

„ . , dn ( 1 \ sin OdOdip , , 

r WsF = ( 1 +3°-0 _ sr ie - (2) 

The corresponding directional diagram is depicted in Fig. QJi. The asymmetry of angular distribution 
([2]) allows to reconstruct the spin state of muon. Moreover, observing the time evolution T(n, t) (mea- 
suring /iSR histograms) allows to trace the time evolution of the muon spin polarization which depends 
sensitively on the magnetic environment. In this paper, we relate the ^SR histogram T(n, t) and a spin 
tomogram w(m,n,t). 

The other aim of this paper is to consider the spin system of muonium Mu = /i + e~ from the tomo- 
graphic point of view. Muonium is formed in certain materials (for instance, in semiconductors) when 
a positive muon /x + picks up an electron. If this is the case, muonium turns out to be a light isotope of 
hydrogen and is of interest not only in solid state physics but also in chemistry [33]. Thanks to a hyperfine 
coupling, the information contained in the electron spin can be passes on to the muon spin which is then 
analyzed via different //SR techniques. Employing the tomographic-probability representation, we are 
interested not only in reconstructing density matrix of muon-electron spin state but also in formation, 
controlling, and detecting entanglement of this system. Such an approach would be useful to realize in 
practice a challenging idea: check Bell inequalities with the help of true spins states not polarizations of 
photons. 

The paper is organized as follows. 

In Sec. [21 we review briefly the tomography of a single spin and derive the relation between the /iSR 
histogram and the qubit tomogram. In Sec. El we focus our attention on muonium and muonium-like 
systems and address the problem of two-spin tomography. In Sec. HI using the tomographic-probability 
representation, the entanglement phenomenon, Bell-like inequality, and positive partial transpose are 
discussed. In Sec. a time evolution of tomograms is considered. In Sec. El we address muonium and 
muonium-like systems in order to introduce their spin Hamiltonian operator. In Sec. [71 the time evolution 
and entanglement of Mu and Mu-like systems are investigated in the tomographic representation. Also, 
a method for reconstructing muon-electron spin state is discussed. In Sec. El conclusions are presented. 



2 Muon Spin Tomography 

Any spin-j state given by the density operator p can be alternatively described by the following 
unitary spin tomogram w^'{m,u) [27j : 

w^'(m, u) = (jm\upu*\jm), (3) 

where \jm) is an eigenvector of angular momentum operators J z and J 2 , m is the spin projection on in- 
direction (m = — j, — j + 1, . . . , j), u is a general unitary transformation which is given by (2j + l) x (2j + l) 
matrix u E SU(N), N = 2j + 1 in the basis of states \jm). The inverse mapping w^'(rn, u) is also devel- 
oped in [27] and implies integration over all unitary matrices u G SU(N) via a corresponding measure. 
We emphasize the unique properties of the unitary spin tomogram ([3]): non-negativity w^'(m,u) > 
and normalization Ylm=-j w^'(m,u) = 1. The physical meaning of the tomogram ([3]) is the probability 
to obtain the spin projection m after a unitary rotation upu^ of the spin-j state is fulfilled. 

For our purposes, it is important to consider the case u £ SU(2), when the unitary spin tomogram 
transforms into a spin tomogram w(m, n) introduced in [23, 24J . Indeed, if this is the case, it is possible 
to parameterize u = R(n) by unit vector n = (cos if sin 9, sin ip sin 9, cos 9) as follows: 

R(n) = e- i(n± - 5)e , n ± = (- sin 0, cos 0, 0). (4) 

For example, in case of qubits (j = 1/2) the matrix of operator @ in the basis of states \jm) is 

R(n) = ( C ° SW2 > --W2)e-<M (5) 

V ; y sin(0/2)e^ cos(0/2) J V ; 

Thus, we obtain the spin tomogram w(m, n) = (jm\R(n)pw (n)| jm) which has a sense of probability to 
obtain m as a result of measuring spin projection on direction n. The inverse mapping of spin tomogram 
w^'(m, n) onto the density operator p reads 

■ 2tt it 

p= Y^ J- [ di P [ sin 9 d9 w( j \m,n(9, if)) D( j) (m,n(9,(p)). (6) 

m=-j 

where D^'(m,n(9,ip)) is a quantizer operator. According to [23,24J, the quantizer operator is expressed 
through Wigner .D-function D^J (a, /3,j) as follows: 

b®(m,n(9,<P)) = (-l)^E E ( 2 i3 + l) 2 E (-l) mi ifi(^,7) 

J3=0m3=- j3 m^m'j.m^-J 

v / j 3 33 W J J J3 V. ,w. / ,-x 

x . / / \jm 1 ){3m 2 \. (7) 

V mi —mi u / y mj —"1-2 m 3 / 

We cannot help mentioning that other different forms of the quantizer operator are also known [25|J34|. 

2.1 Tomography of muon spin state 

As the tomographic technique is to be applied to muon spin, let us consider the special case j = 1/2. 
To avoid confusion, we will denote muon spin jw = 1/2 rather than j, muon spin projection rrv^' = ±1/2, 

and angular momentum operators of muon spin are designated by 3^' = (Jx , J% !> Jz ) = nd", where 
<x = (a x , a y ,a z ) is the set of Pauli operators. 
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Figure 1: The relation between the angular distribution r(n(#, ip)) of positrons e + from muon /i + decay (a) and 
the muon spin tomogram w^(m^ — +1/2, n((9, tp)) of the state | ft) = |j^ = 1/2,77-1^' = 1/2) (schematically 
depicted by an arrow). 



Using definition ^ and formula ©, it is not hard to see that the spin tomogram of a single muon 
reads 

w (m) ( m (A0 , n ) = - + m M Tr [p^ (n • 6")] , (8) 

whereas the quantizer operator can be written in the following form |25j : 



-J + 3m M (n-o-), 



(9) 



with I being the identity operator. 

In order to compare the muon tomogram w^> (m^> , n) and the angular distribution of decay positrons 
r(n), let us consider the tomogram of the state | ft) = \j^' = 1/2, rrv^' = +1/2). Indeed, in this case we 
obtain 



(m) 




u4 M) (+l/2,n(^)) = -(l + cos#), 



w 



(**)< 



1 



■l/2,n(e,<p)) = ~(l-cos9) 



(10) 



The tomogram u>i (+1/2, n(6>, ip)) is depicted in Fig. [lb. Comparing ([2]) and (fTOj) . we derive the following 
relation between the angular distribution of decay positrons (averaged over all positron energies) and the 
muon tomogram: 



„,(M)( m 00 = +1/2, n) = -r(n) - 1, 



u,(M)( m M =_l/2,n) 



2- -r(n). 
2 V ; 



(11) 



The relations (|lip are valid for an arbitrary moment of time t and they give a direct way of transforming 
/iSR histograms into muon spin tomograms w^'(m^> ,n,t). Moreover, it can be easily checked that the 
relations (jlip are adequate not only for pure states (muon polarization |P| = 1) but also for an arbitrary 
mixed state (|P| < 1). As far as negative muons [i~ are concerned, the relation between tomograms and 
angular distributions of decay electrons is readily obtained by exchanging right-hand sides of Eqs. (jlip . 
In the real experiment, we should take into account the spectral responsivity of positron detectors, and 
this might result in changing coefficients of relation (11 If) but linearity of the relation is to be conserved. In 



particular, if there were available detectors sensitive only to maximum positron energy E mSiX = 52.83 MeV, 
it would be possible to measure muon tomogram directly because in this case w^'irw^' = +1/2, n) = 
r Smax (n) and t»( m to = -1/2, n) = 1 -T En ^(n). 

It is known that the muon 2x2 density matrix depends on three real parameters only (vector of 
polarization P) and can be written in the form p = ^(1 + (P • <r)). The muon tomogram w^>(m^' , n) 
is a function on sphere nG5 2 (and a discrete spin projection). Thus, the muon tomogram contains full 
and even redundant information about the system. Nevertheless, reconstruction of the density matrix by 
integration over sphere (see formula ([6])) exploits as much experimental information as possible. However, 
we admit that the practical realization of such an experimental setup (when detectors cover 4tt steradian) 
is a challenge. To eliminate this difficulty one can resort to spin tomography with a finite number of 
rotations [28]. In case of muons, such a tomography exploits only 3 probabilities w^'{m^' = +1/2, n^), 
k = 1,2,3. As it is shown in [28], the inverse mapping onto the density operator exists whenever triple 
product (ni • [112 x 113]) 7^ and reads 



U 3 



I + J2 [2w (M) (m M = +1/2, n fc ) - l] (<r ■ l fe ), 



P= 2 

fe=l 



:i2) 



where vectors 1& form a dual basis with respect to directions {n/ c }^ =1 : 

[n 2 x n 3 ] [n 3 x m] [n x x n 2 ] 



(ni • [n 2 x n 3 ]) ' (ni • [n 2 x n 3 ]) ' " (ni • [n 2 x n 3 ]) 



(13) 

It is worth noting that, from experimental point of view, the advantage of formula (112f) is that it does 
not restrict us to using orthogonal directions only. Though if the errors are presented in measured 
probabilities w^'{rrv^' = +1/2, n^), k = 1,2,3, the reconstructed density matrix is shown to be the 
least erroneous when directions {nfc}^ =1 are mutually orthogonal [28]. In up-to-date experiments the 
background fraction of detector clicks is usually less that 1% of the initial number of counts. At longer 
timescale t 3> r^ the background noise does not allow to measure the muon tomogram because the 
number of muon decays becomes exponentially small. 

3 Muonium and Two-Spin Tomography 

Dealing with muons, we should not forget their physical implementation and behavior in matter as 
well as a concrete experimental setup for muon decay investigation. As muons are to be stopped for 
measuring angular distribution of decay particles, it is important to study behavior of muons in solid 
state materials. On the other hand, the study of /zSR histograms (muon spin tomograms) gives us a useful 
information about the matter in question. Special attention is paid to the investigation of nanostructures, 
and this is why we do not consider muons in gases and restrict ourselves to solid state materials. 

3.1 Muons in Matter 

If negative muon \x~ is incorporated in material, it behaves as a heavy electron and is readily captured 
by atoms to Is ground state. During its way 'down', muon usually loses the major part of its polarization. 
In Is state, a space part of the muon wave function is nonzero at an atomic nucleus. This fact may lead 
to a nuclear capture of the negative muon. As a result, negative muons in matter exhibit relatively 
small degree of polarization and much shorter lifetime (substantially less than r«). It is the reason why 
negative muons (^"SR technique) are not as often used while investigating local matter properties as 
positive muons. However, if a negative muon is in Is ground state inside an atom with zero-spin nucleus, 
the whole system "atom+//~" can be regarded as an acceptor center (see, e.g. [35]). The pseudonucleus 
"nucleus+/i~" has spin 1/2 and keeps a part of muons initial polarization. For example, negative muons 
in Si can form ^Al acceptor centers. The effective moment of the electron shell in such centers can take 



values greater than 1/2, for example, the case j = 3/2 is discussed in the papers [36],[37] . Such systems 
are also of interest for us and we will consider them in Sees. 16.21 and 171 

In contrast, in a condensed matter, positive muons p + are repelled by atomic nuclei and are retarded 
gradually down to Bohr velocities, and then to thermal velocities within the total time interval t ~ 1CP 10 s. 
The retarded muon can either come to rest at an interstitial site of high symmetry between the lattice 
ions (diamagnetic state) or pick up an electron (formation of muonium Mu = p + e~). A spur model of 
muonium formation is discussed in [38]. In addition, muon can also become a part of a chemical bound, 
for instance, with oxygen 0~/i + similarly to OH bound. We will focus our attention to muonium. It 
is worth noting that muonium has almost the same Bohr radius and ionization potential as hydrogen. 
Also, utilizing ultra-low energy muons (with energies from 0.5 to 30 keV) and selecting their velocities, 
it is possible to investigate muoniums formed at the implantation depths from less than a nanometer to 
several hundred nanometers [39[|40], This fact is crucial for studying surfaces and thin films [41J.I42J. 

Muonium is of special interest for us because it can be considered as a two-spin system. Indeed, 
averaging over space part of the common wave function, muon spin (j^ = 1/2) and electron spin (j( e > — 
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1/2) are coupled by hyperfine interaction H^f 
where uj = 4453 MHz. 

The aim of the subsequent subsection is to consider two-spin tomography and apply it to the spin state 
of muonium and muonium-like systems (we will refer to the muon spin 1/2 together with the arbitrary 
effective electron spin p e > as Mu-like system). 



3.2 Two- Spin Tomography 

Let us consider a system comprising two spins: muon spin jW and an effective electron spin j( e > . 
The state of the system is described by a N x N density operator p, where N = (2j^' + l)(2j( e > + 1). If 
we introduce a basis of system states {|*)}i=i, the density operator can be written in the form of density 
matrix pij = (i\p\j). There are two common ways to choose basis states: 

(i) eigenvectors \jw rn^')\p e ' m^ e ') of individual angular momentum operator J z <g> J z , 
(ii) eigenvectors \LM) of the total angular momentum operators L z and L, where the total spin L 
can take values L = \j^> — j^\, . . . ,j^' + j^ . 

Evidently, these two bases are related by the matrix Ucg composed of the Clebsch-Gordon coefficients 

C L ( f )mM . (e)m(e) = (j^) m Mj(e) m (e)\LM). For muonium (j^) = 1/2, j^ = 1/2) the matrix U C g reads 
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(14) 



Consequently, in the probability representation muonium can be described by two tomograms, namely, 
an individual two-spin tomogram w{rn^\r\S 11 ' ,rn^ e 'rv- e >) and a total two-spin tomogram u/ >(M, N) 
discussed below. 



3.2.1 Individual Two-Spin Tomogram 

Extending the notion of spin tomogram ([3]) to the case of bipartite spin system, we immediately 
obtain the individual unitary two-spin tomogram defined through 



w(m^\m^,U) = {jV>mM,j®mM\U p U^\j M m^\j {e) m^) 



(15) 



where U is a general N x N unitary transformation, A^ = (2j^> + l)(2j( e ' + 1). If U is factorized into 
U = -R(n^) <S> R(n.( e >), where the rotation operator R is given by formula (j4]), we readily obtain the 



so-called individual two-spin tomogram of the form 

w(m M ,n W ,m (e U (e) ) = (j M m^ ,j^m^\R(n^) (S) R{n^) p R^n^)® R^ n V)\jM m 0*)j<fi) m {*)) } 

(16) 
which has a clear physical meaning: the individual two-spin tomogram is the joint probability to obtain 
the value m Ji ' in measuring the muon spin projection on direction n^ and the value m^ e > in measuring 
the electron spin projection on direction r\> e >. Here, measurements are assumed to be simultaneous. 
Though, in this case, such measurements are not prohibited by quantum mechanics, it is not evident how 
to realize them in practice. The reconstruction of the density matrix of the whole system is 

j(p) , . j( e ) , , 

p= J2 f^- £ j^r w(^ M ,nM >TO W jn W ) o&W) (m M in M )s ^CiW) (m W in (e) )| 

(17) 
where the local quantizer operator D^'(m,n) is given by formula (J7J). 

As mentioned above, the practical realization of measuring the individual spin-tomogram is a chal- 
lenge. In the fact, the //SR-method (discussed in Sec. [2]) allows measuring the muon part of the tomo- 
gram only. Using the probability interpretation of the tomogram (116p . we can say that only the following 
marginal distribution is measurable via pSH technique: 

j(e) 

w(m^\n^) = J2 w(m (M \n^^m (e) n( e )) = Tr[p#(nW)|jWm ( ^)0' M m ( ^|^(nW)(8)jW] 

m( e )=-j( e ) 

= {j' M m W |fl(n W ) Tr e [p] R\n^)\j^m^), (18) 

which is nothing else but the muon spin tomogram of the reduced density operator Tr e [/5] (trace is taken 
over electron spin states). Here, we used the following resolution of identity: 

/( e ) = J2 RHn {e) )\J (e) ™ {e) )(J {e) ™ {e) \R(n {e) ), (19) 

m(e)=-j(e) 

which is valid for any direction m e ). Since the summation over spin projections m^ e > leads to exclusion 
of \v- e > in the reduced tomogram (I18p (this phenomenon is also known as non-signalling), we draw a 
conclusion that it is impossible to reconstruct the muonium density matrix by measuring the reduced 
tomogram (|18p only. 

However, if we consider the individual unitary two-spin tomogram (fT5j) . we will see that the reduced 
tomogram 

w{m^\lJ)= Yl w(m itl \m ie \lf) = {j^ ) m M \Tr e [U ptP]\jM m M) (20) 

m (e)=-j(.e) 

does contain the full information about the muonium spin state. The physical meaning of the distribution 
function fj20f> is a probability to obtain the muon spin projection on z-axis being equal to m^> after a 
unitary rotation U is fulfilled over the common density operator p. In Sec. 17. 1[ we will outline how to 
reconstruct the muonium (and Mu-like) density matrix by using measurable muon tomogram (I20p . where 
the role of unitary transformation U is played by an evolution operator ^(t). 

3.2.2 Total Two- Spin Tomogram 

From the mathematical point of view, the total unitary two-spin tomogram w^ L '(M,U) is merely a 
set of diagonal elements of the transformed density operator in the alternative basis of states \LM): 

w iL) {M,U) = (LM\U pU ] \LM). (21) 



Nonetheless, the tomogram w { '■ '(M,U) has a physical meaning different from that of the individual 
unitary two-spin tomogram f)15|) . Actually, w^ L '(M,U) is the probability to find the system in L-spin 
configuration and to obtain the total spin projection onto z-direction equal to M after the density operator 
p experiences a unitary transformation U. Obviously, the tomogram w^ L '(M,U) is normalized, i.e. 
Sl ^m=-l w (M, U) = 1. The total unitary two-spin tomogram (|2ip contains the same information 
as the individual unitary two-spin tomogram (|15|) because both of them are nothing else but diagonal 
representations of the rotated density operator U ' pW , with the only difference being in the basis used: 
either {|/^m>))|/ e W e ))} or {\LM)}. 

Similarly to the case of individual two-spin tomography, where we specified the unitary transformation 
in the factorized form of individual rotations, i.e. U- m d = R(n^>) (g) R(n^ e '), it is tempting to introduce 
separate rotations of subspaces {\L'M')} and {\L"M")} with L' ^ L" . For this reason, we introduce the 
following operator: 

j(M)_|_j(e) 

L>tot=£ (L,) (n L 0©£ (L "W)®---^ R (L) (n L ), (22) 

where by sign © we designate the direct sum. This means that the matrix representation Utot of the 
operator C7tot in the basis of states {|LM)} is block-diagonal, with the L-th block being the (2L + 1) x 
(2L + 1) matrix # (L) (n L ) defined by formula flU). 

It is of vital importance to trace how the problem of addition of angular moments is related with 
individual and total two-spin tomograms. The result is 

& L \N) = ^(N) © R( L "\N) © . . . = flCJ w )(N) g, j?0' (e) )(N) (23) 

L= |j(M)-j(e)| 

or in the matrix form 

j(M)+j( e ) 

R( L \N) = R( l '\N) © i?( L ")(N) © . . . = C/J G ^' (M) )(N) © i^' (e) )(N) U CG , (24) 

L=\jto-j(°)\ 

where Ucg is a matrix of Clebsch-Gordon coefficients. 

Arguing in the same manner as in case of individual two-spin tomogram f 1 1 6 [) . we introduce the 
following total two-spin probability distribution function of the form 

j(f)_|_j(e) j(f)_|_j(e) 

f( L \M,N) = {LM\ R iL \N) p R.W(N)\LM). (25) 

L=\jM-j(*)\ L=\jM-jM\ 



It turns out that the knowledge of the probability distribution function (125ft is not enough, in general, to 
reconstruct the system density operator of an arbitrary mixed state [33] . Consequently, we cannot regard 
the function f( L >(M, N) as a true tomogram. However, if the information in (j25j) is supplemented by 
probabilities with which pure states appear in the mixed one, the mapping p — > /' \M, N) is proved 
to be invertible [33]. 

Despite the probability distribution (i25j) is not informationally complete, it can still be treated as a 
tomogram of identical spin-j particles |44j . Indeed, in case of a bipartite system with indistinguishable 
constituents, the density operator of the whole system p is known to be invariant with respect to a 
particle permutation. Hence, the density matrix in the basis of states {\LM)} is block-diagonal and can 
be reconstructed as follows: 

P indistinguishable = P^ = E J^ f {L) (M,N)D^ (M,N), (26) 

L=0 L=0 M=-L g 2 



where the quantizer £K '(M, N) is determined by formula ([7]). As far as other aspects of the addition 
of angular moments in the tomographic-probability representation are beyond the MuSR-method, they 
will be discussed elsewhere. 

In order to conclude this consideration, we derive the relation between the individual two-spin tomo- 
gram (|16p and the total unitary two-spin tomogram f)21j) 



w 



(L) 



(M, U)= Y, J 



(m)=_j(m) a2 



dnM ' r cW e ) 

Air 2-^ J Air 

m( e )=— j( e ) g2 



u)(m w ,n M ,m (e) ,nW) 



x(LM\U D (jW \m^\n^)®D( j{e> \m( e \n^) tf\LM). 



(27) 



4 Entanglement and Separability in Probability Representation 

A state of two-spin system (muon and electron spins) is called separable if the density operator p of 
the whole system can be resolved into the following convex sum 



/^separable / Pi P 



(M) o, -(e) 



o<Pi<i, 5Z« = 1 ' 



(28) 



otherwise the state is called entangled. Since the individual two-spin tomogram io(m'^, n'^^m' 6 ^, n' e ^) 
contains the same information about the system involved, the definition of separability (|28p is readily 
reformulated in terms of tomograms. In fact, the state is separable when 

^se P arabie(mW,nW,m( e ),n( e )) = £> u,^ (m^ nM) W f <e)) (m (e) , n( e )). (29) 

i 

If the resolution (I29p is not fulfilled for any convex parameters pi, then the state is entangled. 

4.1 Entanglement Detection via Bell- like Inequality Violation 

To detect the entanglement of a two-qubit system it is convenient to check the violation of the following 
Bell-like inequality (see, e.g. |45j): 

\B\ < 2, 



(30) 



where the Bell-like number B reads 



B = Tt 



f w(+, n^,+,nS e) ) w(+,n^,+,n { 2 e) ) 



J 
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(eh 



w(+,n^,-,n?>) w(+,n^ 



» 



» 



,(e) 



w{-,n^,+,n?>) w(-,n^',+,n 



(m) 



(/■') 



\ w(-,ii y f',-,n 



.(e), 



?(' 



» 



4 e) ; 






ir( 



V) 



, n 2 , -r, ±± 1 



+,n 



(eh 



w(-,n£ 



(m) 



.(<0> 



^(+,n^,+,4 e) ) \ 



(31) 



where the sign ± denotes spin projection ±1/2, respectively, and the matrix J? is defined through 

/ 1 -1-1 1 \ 

1 -1-1 1 

1 -1-1 1 

-11 1-1 



J 



\ 



(32) 



I 



If state is separable, then the inequality (|30p is necessarily valid. Thus, its violation is a direct evidence 
of the system entanglement. Note, that we do not resort to the density matrix formalism while dealing 
with such an approach. 



4.2 Positive Partial Transpose in the Tomographic-Probability Representation 

Two-qubit and qubit-qutrit states are known to be separable iff positive partial transpose p ppt of the 
system density matrix p results in a new density matrix (Peres-Horodecki criterion [46,47j). If p ppt is not 
a density matrix, this is the direct evidence of the entanglement. The Peres-Horodecki criterion differs 
substantially from the approach outlined in the previous subsection because it provides the necessary 
and sufficient condition for a two-qubit state or qubit-qutrit state to be entangled whereas the violation 
of the Bell-like inequality (130|) is only a sufficient but not necessary indicator of entanglement. 

Positive partial transpose (PPT) of bipartite state p leaves the second subsystem undisturbed and 
makes a transpose of the first subsystem, which is equivalent to the complex conjugation of the first 
subsystem (time reversion). Since the complex conjugation of qubit density matrix p = ^(I + (P • <x)) 
boils down to the replacement P y — > —P y , the corresponding qubit tomogram w(m,n) transforms into 
w(m, n ppt ), where n ppt = (n x ,—n y ,n z ). Consequently, in case of two qubits (muonium) and qubit- 
qutrit system (muonium-like system) , the action of positive partial transpose on the individual two-spin 
tomogram (|16p is 

w ppt (m^ , nW , m^ , n^ ) = w(m^ , n^ ppt , m^ , n( e ) ) . (33) 

In the density- matrix formalism, the following step is to check a positivity of the matrix p ppt . It is shown 
in [38] that eigenvalues {A,}^ =1 of a 4 x 4 Hermitian matrix A with unit trace (TrA = 1) are non-negative 
iff M<i = Yl AjAj > 0, M3 = Yl AjAjAfc > 0, and M\ = ^ XiXjX^Xi > 0. Using the trace equality 

i<j i<j<k i<j<k<l 

J2 A« = 1, these quantities can be expressed through traces TrA n as follows: 

i 

M 2 = i (1 - TrA 2 ) > 0, (34) 

M 3 = - (1 - 3TrA 2 + 2TrA 3 ) > 0, (35) 

— (\ -fiTrA 2 +3rTrA 212 
24 

Substituting p ppt for A in (f34"|) - (f36|) . it is possible to detect unambiguously the entanglement of two-qubit 
system by observing a violation of at least one inequality. The value of M2 turns out to be invariant with 
respect to PPT transformation. In view of this, we introduce the following entanglement measure E for 
two qubits: 

E = \M 3 \ + \M 4 \ - M 3 - M 4 . (37) 

In the tomographic-probability approach, the values of M3 and M4 read 



M 4 = — (l - 6TrA 2 + 3 [TrA 2 ] 2 + 8TVA 3 - 6TrA 4 ) > 0. (36) 



m 3 = i Yl J2 K pt - 3u,ppt * wPPt + 2wPPt * wPPt * ^ ppt ] ( m(M) > n(M) > m(e) > n(e) ) > ( 38 ) 

M 4 = -U 3( Yl ^K Pt « Ppt ](™ ( "U W ,m (e U (e) )) 2 + Yl Y^ wPPt - 6wPPt * wPPt 

+8w ppt * w ppt * w ppt - 6w ppt * w ppt * w ppt * w ppt ] (m M , n^ , m (e) , n^ ) j , (39) 

where by * we denote a star product with the kernel 

2 
K(m in]m2 n 2) wiinjmjM^, rni~n-im2Tn.2) = I I 7 + 9m^77i" (n^ • n'/) + 3mjm,j (n 4 - • nQ 

4 = 1 

+3m 4 mf (rij • n") + il&mim'jm'l (nj • [n^ x n'/]) . (40) 



Thus, the measure (J37J1 is expressed through spin tomograms only. 
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5 Evolution of Tomograms 

Evolution of pure states l^) is governed by Schrodinger equation ih-fijf- = H\ip), where H is a 
Hamiltonian of the system. If we consider generally mixed state with the density operator p, the evolution 
obeys the von Neumann equation ih-S = [H, p] which is similar to Liuville equation in classical statistical 
mechanics. Taking advantage of mapping p(t) — > w(m,u,t) for a single spin-j particle (Eq. fl3J)), the 
von Neumann equation with the initial condition takes the form 

OW ( TYl U t] 

ifh ' ' =[fH*w-w*f H ](m,u,t), w(m,u,t = 0) = w (m,u), (41) 

where fjj (m, u, t) is a tomographic symbol of operator H obtained by replacing p — >• H in formula (|3J); a 
star * implies the so-called star product of tomographic symbols. In other words, 

[fa * w](m,u,t) = 2_. I du' I du" /#(m' ,u ,i)/jy(m / ,u',t)K(m',u',m",u" ,m,u), (42) 

am' wi" i 



m ,m' 



where J du is an integration over the corresponding measure, K(m' ,u' ,rn",u" ,m,u) is the star-product 
kernel of the form K = Tr[D(m' ,u')D{m'',u'')U(rn,u)\, where D and U are quantizer and dequantizer 
operators, respectively. Explicit form of the spin tomographic star-product kernel is found in [25.26JI34"]. 
Formal solutions of Schrodinger and von Neumann equations are expressed through the unitary evo- 
lution operator &(t) as follows: 

\m) = #(*)W0)>, P(t) = &(t)p(0)&Ht), (43) 



W(t) = Texp 



■lll k[T)d \ 



(44) 



where T is a time-ordering operator. 

Using the relation (|43p as well as the mapping ([3]), it is not hard to see that the time evolution of a 
unitary spin tomogram reads 

w(m, u,t) = wo(m,u^(t)). (45) 



which is nothing else but the formal solution of the tomographic evolution equation (|4ip . 

As far as muonium and muonium-like systems are concerned, the evolution of the individual unitary 
two-spin tomogram (|15p reads 



w(m M ,m {e \U,t) = wo(m M ,m {e \U&{t)), (46) 

whereas the evolution of the individual two-spin tomogram (|16p has a rather complicated form 
u;( m 0*),n^,m (e) ,nW,t) 

J r Jfj(W { I" Wn( e ) 

m(M)=-j(M) 5 2 m( e )=-j( e ) 5 2 

xt(t)^' ( " ) )(mW,nM) ® ^''^(mfV^f *(t) flt( n W) 8 #(nW)|jM m M, jW m W). (47) 
6 Muonium and Mu-like Systems: Hamiltonian 

Hamiltonian of a muonium-like system has the following most general form: 

^general (ff, s) = flbfa) + V(q, s), (48) 
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where q and s are responsible for position and spin variables, respectively; Ho(q) is a part depending 
on space coordinates only; and V(q, s) is a part containing spin variables. For muonium-like systems, 
Coulomb and spin-orbital interactions are much greater than the hyperfine interaction {Hq 3> V). For 
this reason, in our case the operator V(q, s) will only contain the hyperfine interaction and an interaction 
of magnetic moments Ha and fj, e with local magnetic fields. On averaging over space coordinates, the 
operator {ip(q)\V(q, s)\ip(q)) acts on spin degrees of freedom, i.e., it can be represented by N x N spin- 
Hamiltonian operator H with N = (2j^ + l)(2j^ + 1). 



6.1 Muonium Hamiltonian 

We will consider three main kinds of muonium Hamiltonian: 

• Free muonium atom in vacuum is governed by hyperfine interaction 

H M = hcj (J^ J (e) ) = ^o(4 M) 4 e) + jf ® 4 e) + J^ Ji e) ), 



(49) 



where the frequency of hyperfine interaction ojq = 4453 MHz. It is worth mentioning that such a 
Hamiltonian provides an adequate description of muonium not only in vacuum but also in quartz 
and ice. Frequency of hyperfine splitting in such materials is almost the same as in vacuum within 
the accuracy of the experiment. The correction caused by interaction of the electron quadrupole 
electrical moment with inhomogeneous crystal field can lead to anisotropic Hamiltonian (in quartz 
at temperatures T < 100 K). 

Isotropic Hamiltonian of a muonium-like system in solids in presence of magnetic field B 

H M u = A(jM JW) - g^(B ■ jM) / (e) + geHel^ ® (B • j( e )), (50) 

where in general A ^ Ttujq, g^ = 2, [i e is the Bohr magneton, and (B • J) = B X J X + B y J y + B Z J Z . 

Anisotropic Hamiltonian of anomalous muonium (Mu*) which is present, for instance, in Si and 
many other semiconductors with a diamond structure of crystal lattice 



H Mu * = H Mu + AA(m • j(">) (01 • j (e) ; 
where Ot is a unit vector determining the axis of axial asymmetry. 



(51) 



If the Hamiltonian is time- independent (magnetic field is steady), then the corresponding unitary 
evolution operators can be found analytically. In the basis of states \j^'m^> ,j^ e >m^ e '), the unitary 
evolution matrix of hyperfine interaction reads 



(2 



^hf(t) 



± — iwo*/4 











o\ 



1 + e iuJot 1 - e iuJot 
1 - e iuJot 1 + e iuJot 
2 



(52) 



In order to find the evolution matrix of Mu in presence of magnetic field, the intensity B is chosen to be 
aligned along z-direction, i.e. B || z (longitudinal field). In this case, we obtain 



/ -i(2a-b-)t 



%lu z(t) 



Aat 











cos ct + i (6+ /c) sin ct —i(2a/c) sin ct 

—i(2a/c) sin ct cos ct — i{b + /c) sin ct 

e" 






\ 












i(2a+b- 


"J 



(53) 
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where a = A/ (Ah), b± = B(g^^ ± g e fi e )/(2H), and c = yj 'A 2 + B 2 (g fJi iJ, fJi + g e ^ e ) 2 /(2h). 

If the magnetic field is aligned along x-axis (transversal field), the unitary evolution matrix ^Mu x(t) 
is symmetrical, and its matrix elements read 



[%Aux(t) 
[%u x(t) 
[%u *(<) 
[%u,(t) 



11 



22 



12 



13 



[%u x(t)] M 
[%u *(*)]33 
[%u x(*)]34 
[^Mu x(t)]24 



[e" iaf cos 6_t + e iat (cos ct - z(2o/c) sin ct) 
[e" ia * cos 6_t + e ia< (cos ct + z(2o/c) sin ct) 



[e- iat sin b-t + (6+/c)e mt sin ct] , 



- [e- iat smb_t - (b + /c)e mt sin ct] , 



1 



u = 2 [e 



—iat 



cos 6_t - e m< (cos ct - i(2a/c) sin ct)] , 



1 



23 



[e" iat cos 6_i - e ia< (cos ct + i(2a/c) sin ct) 



(54) 
(55) 
(56) 
(57) 
(58) 
(59) 



Finally, if direction of the magnetic field coincides with y-axis, the unitary evolution matrix % u y (t) is 
expressed through matrix % u x (t) as follows: 



%u y(t) 



I [%u,(t)]ll -i[%u *(*)]l2 -i[^Mu x(*)]l3 -[%u,W]l4 \ 

J[*Mux(i)]21 [%ux(t)] 2 2 [^Muz(i)]23 -«[%o ,W] 2 4 

i[%u *(t)] 3 l [%ux(t)] 3 2 [%xW]33 ~i[%u x(*)]34 

V -[ < 2fo«*(t)]41 »[% os(t)]42 i[%u,(f)]43 [%u,W]44 / 



(60) 



As concerns the anisotropic hyperfine interaction (|5ip . we will present below the explicit analytical 
solutions for three cases. The first one is OX \\ z and B || z, in this case we obtain the following unitary 
evolution matrix ^Mu*zz(t)'- 



°fa 



Mu*zz 



(t) 



J,(a+d)t 



/ e -i(2(a+d)-b-)t 
















cos ct + i (6+ /c) sin ct —i(2a/c) sin ct 

—i(2a/c) sin ct cos ct — i(b + /c) sin ct 

e -i(2(a+d)+b-)t 



\ 



where d = AA/(4h). 

In case OX \\ x and B || z, we obtain symmetrical matrix % u » K (t) with the elements 

[^Mu^(i)]n=e- iat [cos/t + i(6_//)sin/t], 

[&Mu*xz(t)]l2 = [%An*xz(t)]t3 = [&Mu* xz(t)U = [%u*, 2 (t)] 3 4 = 0, 

[%An*xz(t)]u = -ie~ iat (d/f) sin ft, 
[%Au*xz(t)}22 = e iat [ cos ht + i(b+/h) sin ht] , 
[%Au* xz (t)}23 = -ie iat [{2a + d)/h] sin /it, 
[%An*xz(t)]aa = e iat [ cos ht - i(b+/h) sin ht] , 
[%An*xz(t)]u = e~ iat [ cos ft + »(&_//) sin /t] , 



(61) 



(62) 
(63) 
(64) 
(65) 
(66) 
(67) 
(68) 



where / = y/{AA) 2 + 4B 2 (g^ - g e ^e) 2 /(^h) and fc = ^/p + AA/2) 2 + 5 2 ( 5/ ^ + g ef i e ) 2 /(2h). 

Finally, in case OX \\ y and B || z, we obtain the symmetrical unitary evolution matrix ^Mu*yz(t) whose 
elements coincide with elements of the matrix % u * IZ (t) except for [^Mu*j/z(*)]i4 = ~~ [^Mu*a;z(^)]l4- 
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6.2 Mu-like Hamiltonian 



The muonium-like systems are described by Hamiltonians ()49[) — f)5ip . where the effective moment of 
the electron shell j( e > is greater than 1/2, i.e. the dimension of the operator j( e > is greater than 2. In the 
particular case j( e > = 1, the effective hyperfine interaction will result in the following unitary evolution 
matrix: 



%& 



-like 



(t) 



( fx 
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y + 
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n) 



(69) 



where matrix elements "V are expressed through the strength of hyperfine interaction A 






-iAt/{2h) 



%, 



-ie 



iAt/(AK) 



2V2 / 3) sin(3 At /(4h)), 



JAt/(4h) 



[cos{3At/{4h)) ± (i/3) sin(3i4t/(4R))] . 



(70) 

(71) 



The entanglement of a Mu-like systems can be detected by using tomogram w(mM 1 ' , lay 1 ' ,rnS e > , n' e )) 
together with the qubit portrait method |45j and Bell-like inequality (|30p . Moreover, PPT inequalities 
(|34 p ~ (|36p are necessary conditions for separability, so the non-zero value of quantity (|37p is a direct 
indicator of the entanglement of a Mu-like system as well. 

It is worth noting, that we restrict ourselves to the unitary evolution of muonium-like systems and do 
not consider relaxation processes. Their influence on system behavior and entanglement will be discussed 
elsewhere. 



7 Evolution and Entanglement of Mu-like Systems 

In this section, we exploit the developed tomographic approach to solve the problem of evolution 
and entanglement of Mu, Mu*, and Mu-like systems. We are going to indicate entanglement of Mu-like 
systems and propose experimentally accessible techniques for creation of entangled states. These two 
points can be used in order to check experimentally Bell inequalities with solid state spin states. 

Let us consider the muonium atom. 

As it is shown in Sec. El for a given Hamiltonian, the evolution is thoroughly determined by the 
initial state of a system. In order to demonstrate how the developed method works, we focus on the 
simplest and physically clear initial state. Since the retarding time of muons in solids is ~ 1CP 10 s, the 
muons conserve their polarization and are assumed 100% polarized. In contrast, during the formation of 
muonium the captured electrons have either up or down spin-z projection, i.e. the electron subsystem 
is in the maximally mixed state (non-coherent mixture). The reference time t = is the formation of a 
muonium. The initial state of a muonium is 



p(0) = I W I ^/ (e) . 



(72) 



Consequently, the initial individual two-spin tomogram (|16p reads 



wo(m M ,n^,m^,n^) 



1 (\ 



2 V2 



+ m W nW 



(73) 



that is the state is factorized (simply separable) and tomogram (|73p does not depend on the variables m 



(e) 



and rr e ) ascribed to the electron. Substituting the initial tomogram ([73 p for 



wo{m 



00 



riW.mW. 



n 



(e) 



) in 
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Figure 2: Time evolution of the reduced tomogram w(l/2, n^\ t) of Mu in quartz (left-up) and Mu-like system 
with effective electron moment j( e > — 1 (right-up). Curved lines correspond to n^ || z, horizontal (fine solid) lines 
correspond to n^ || x and n^) || y. Time dependence of entanglement for Mu in quartz (left-bottom) and for 
Mu-like system (right-bottom). The initial state is chosen to be {2j^ + 1) _1 | if){it I ® 1^ ■ Longitudinal magnetic 
field B || z with magnitude B = (heavy solid line), B — 790 Gs (dashed line), B = 1580 Gs (chain line), and 
5 = 3160 Gs (dotted line). 



(|47|) and separating variables {m^' ,11^'} and {m( e ' ,n( e >}, we readily obtain a solution in the form of 
the time-dependent tomogram 

1 



w(m ( ^,n^,m (e) ,n( e ),t) 



1 + (jW m W, jW m W\R(nW) R(pS e >) 

X$(t)&M ® I^^{t) £t( n (*0) ^ flt( n (e))|jO*) m (A*) 9 j(e) m (e)\ 



(74) 



which is suitable to describe the evolution of muonium governed by any unitary evolution f)52[> — (|62p . 
Let us first consider the evolution in approximation of hyperfine interaction with Hamiltonian (|49p . 
On substituting the unitary evolution matrix (j52|) for fy (t) in ([74")) , the direct calculation yields 



w{m 



(P) n M ™( e ) n( e ) 



n 



m y 



,nW,t) 



jfl + mWnM +m (e) ni e) + (m M n^ - m (e) 4 e) ) cos w t 
+2m M m (e V (M) x nW], sin wot). (75) 



As mentioned in Sec. 13.2.14 the MuSR-experiment enables us to measure the reduced tomogram (to- 
mogram of muon subsystem). In case of the hyperfine interaction evolution, the reduced tomogram 
reads 
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w{m^\n^\t) = - \l + m^n^\l + cos w t) 



1 

2 L 

The last term in formula (|75p is responsible for 
entanglement of muon and electron spins. The direct 
calculation shows that the maximum Bell-like num- 
ber (j3"0|) for tomogram (fT5j) is max ( M , e ) B = | sinu;o£|. 

n l,2 

Thus, the Bell number satisfies \B\ < 1 and is unap- 
propriate for detecting entanglement of free muonium. 
Partly, it is due to the initial state ([72l) ~ ([73j) which is 
obviously mixed whereas the operation of the Bell- 
like number is best in case of pure states. Applying 
the PPT criterion to the tomogram (|75p . we obtain 
E = (1/128) sin 4 ujQt, which reveals the entanglement 
of a free muonium at the time moments t ^ nk/coQ, 
k = 0, 1, 2, . . . Time evolution of the reduced tomo- 
gram and the entanglement measure is visualized in 
Fig- [2] (muonium in quartz, zero field). 

Dealing with more complicated systems, for the 
sake of brevity, we will omit analytical solutions of the 
tomographic evolution equation (147p and calculations 
of the entanglement measure ()37[) . 

In Fig. [2J we present the reduced (muon) tomo- 
gram and the entanglement measure for the simplest 
Mu-like system "muon + spin 1" which is governed 
by unitary evolution (|69j) . 

As far as muonium in quartz is concerned, we ex- 
ploit the Hamiltonian ([50]) with A = (2irh) ■ 4404 MHz 
and different magnetic fields B = 0,B C /2,B C ,2B C , 
where B c = 1580 Gs is a critical field, i.e. B(g e fi e — 
dii^fj.) = A. In Fig. [21 we illustrate the behavior of 
the observed reduced tomogram w{mS^' , n'^), t) and 
the entanglement measure E when B [| z (longitudi- 
nal field). While increasing the magnetic field inten- 
sity B, the tomographic value w(mS IJ '' = 1/2, ng , i) 
oscillates with increasing frequency and tends to 1 due 
to a preferable muon spin polarization along positive 
z-axis and electron spin polarization along negative 
z-axis (separable state). For the same reason, the en- 
tanglement measure E decreases substantially when 
the strength of magnetic filed is greater than the crit- 
ical value B c . The probability to obtain muon spin 
projection +1/2 along x- or y-axis is exactly 1/2 and 
is not sensitive to the magnitude of magnetic field 
aligned along z-axis. 

The case of muonium in quartz, B || x (transver- 
sal field) is depicted in Fig. [3l Both tomographic val- 
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ues w(l/2,n z ,t) and w(l/2,n]f , t) exhibit a rather 
complicated dynamics, whereas w(l/2,riy , t) = 1/2. 
Again, the essential decrease of entanglement is observed when B > B c 



0.4 0.5 

t*10 9 ,s 

Figure 3: Time evolution of the reduced tomogram 
55(1/2, n^,*) of Mu in quartz with n^ || z (up), 
n^) || x (middle, curved lines), and n^') || y (middle, 
horizontal fine solid line). Time dependence of entan- 
glement for Mu in quartz (bottom) . The initial state is 
I ft) (ft I ® \l^ ■ Transversal magnetic field B || x with 
magnitude B = 790 Gs (dashed line), B = 1580 Gs 
(chain line), and B = 3160 Gs (dotted line). 
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Figure 4: Time evolution of the reduced tomogram w(l/2, n^\ t) of Mu* in Si with n^ || z (up, curved lines), 
nW || x and nS^) || y (up, horizontal fine dotted line). Time dependence of entanglement for Mu* in Si (three 
bottom figures). The initial state is | ff) (ff | ® \l^ ■ Anisotropy direction 9t || x. Longitudinal magnetic field B || z 
with magnitude B — (dotted line), B = 10 Gs (dashed line), B — 33 Gs (heavy solid line), and B — 100 Gs (fine 
solid line). 
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Anomalous muonium behaves as an ordinary muonium if 91 \\ z and B || z, with the only difference 
being in time scale. The case 9t || x and B || z is illustrated in Fig. [H we depict the reduced tomogram and 
entanglement evolution for anomalous muonium in silicon (A = 92.595 MHz, AA = —75.776 MHz [49050] ) 
in presence of different magnetic fields B = 0, 10, 33, and 100 Gs (critical field B c = 33 Gs). A two- 
frequency modulation of the reduced tomogram and a decrease of entanglement are observed if B > B c . 

7.1 Reconstruction of the Initial State 

The strong dependence of the reduced (muon) spin tomogram on time can serve as an indicator of a 
system Hamiltonian and allows to determine numerical values of parameters (^4, AA, B, 91). Given the 
evolution operator ^f(i) and the measured reduced tomogram u>(m^,n^, t), it is possible, in general, 
to reconstruct the initial state /5(0) [51) . Indeed, fifteen independent values {w(l/2, n^ , ti)}, k = x,y,z, 
I = 1,2, ... ,5 allow the reconstruction of fifteen independent real parameters of 4 x 4 density matrix 
p(0). The requirement on time moments {ti}f =1 is to avoid a periodical coincidence. The important 
requirement on unitary evolution ^(t) is to break an internal symmetry of the two-spin state, i.e. an 
anisotropic term like AA in (15ip is to be presented in the system Hamiltonian. It is worth mentioning that 
the idea of using time unitary evolution to measure a symplectic tomography of Bose-Einstein condensate 
was successfully applied in the paper [52j. 

8 Conclusions 

To conclude we point out the main results of our work. 

We applied the tomographic probability approach to the problem of muon spin rotation/relaxation/ 
resonance experiments. The relation between MuSR histogram and muon spin tomogram was estab- 
lished. In view of this, all MuSR experiments can be interpreted as measuring the reduced tomogram of 
a quantum state of two-spin system, namely, muon-electron one. The density matrix of this system was 
mapped onto a joint probability distribution of muon and electron spin projections onto their quantiza- 
tion directions. Such a tomogram was shown to completely describe any muon-electron quantum state. 
The reduced density matrix describing the muon state was mapped bijectively onto a spin-tomographic 
probability distribution of the muon spin projection onto a chosen direction (muon spin tomogram). 

Time evolution of muon-electron system was investigated within the framework of the tomographic 
representation. The entanglement phenomenon between muon and electron spins was shown to appear 
due to different interactions including the hyperfine interaction, the anisotropic hyperfine interaction and 
the interaction of muon and electron magnetic moments with an external magnetic field. We studied the 
behavior of the Bell-like tomographic inequality and showed that, for a discussed initial system state, 
there is no violation of this inequality. The PPT criterion of separability was discussed in the tomographic 
probability representation. The negativity indicator of entanglement was applied to the muon-electron 
system. Using such a tomographic entanglement measure, we studied the influence of different kinds of 
the discussed interactions onto the entanglement phenomenon. The detailed analysis of entanglement 
under the influence of dissipation within the framework of kinetic equation with relaxation terms is a 
problem for further consideration. 
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